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Abstract 

For a Gaussian process X and smooth function /, we consider a Stratonovich integral of f(X), 
defined as the weak limit, if it exists, of a sequence of Riemann sums. We give covariance conditions 
on X such that the sequence converges in law. This gives a change-of-variable formula in law with 
a correction term which is an Ito integral of /"' with respect to a Gaussian martingale independent 
Qh ' of X. The proof uses Malliavin calculus and a central limit theorem from [5]. This formula was 

(-H known for fBm with H = 1/6 |10J . We extend this to a larger class of Gaussian processes. 

-I— > 

1 Introduction 

Let X = {X t , t > 0} be a centered Gaussian process, and let / : R — » R be a C°° function such that 
/ and its derivatives have at most polynomial growth. We define the Stratonovich integral of /' with 



> 

respect to X, denoted, 



f(X s )d°X s , 



o 



as the limit in probability of the trapezoidal Riemann sum, 

1 L"^ 1 r w n 



2 

i=0 

when that limit exists. In a 2005 paper, Gradinaru e£ al. ( 5 ) studied this integral, and identified 
conditions on X under which the Riemann sum converges. In particular, the sum converges for any 
fractional Brownian motion (fBm) {B t ,t > 0} with Hurst parameter H > 1/6, in which case the 
following change-of-variable formula holds: 

f(B t ) = f(0)+f t f'(B s )d°B s . (2) 
Jo 

Subsequently, JXOJ examined the end point case H — 1/6. Here, it was proved that (TT]) converges weakly 
to an Ito-like expansion formula, consisting of a stochastic integral and a correction term in the form 
of an Ito integral of /"'. In this case, we have the weak change-of-variable formula, 

f(B t ) k /(0) + I* f(B s )d°B s J f^(X s )dW s , (3) 
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where W is a Brownian motion, independent of B 7 with variance given by (|5|). 

In this paper, we consider the behavior of (UJ for a more general class of Gaussian processes, which 
are characterized by conditions on the covariance. Convergence follows from a central limit theorem 
first proved by Nourdin and Nualart in [5]. This theorem is based on Malliavin calculus, and applies to 
a sequence of multiple Skorohod integrals. In this paper, we give a set of six covariance conditions on 
the process X, which lead to weak convergence of the form ([3]). These conditions are satisified by fBm 
with H = 1/6. As an application, we have found that the conditions are met for three fBm-derived 
processes, including 

• Bifractional Brownian motion (bBm) with parameters HK = 1/6; 

• 'extended' bBm with K e (1, 2) and HK = 1/6; and 

• sub- fractional Brownian motion with parameter h = 1/3. 

In the prequel to this paper ([6]), we applied the same central limit theorem to a 'midpoint' Riemann 
sum of the form 

f'(X aj-i )(X2j — X ij-i ). 

i-i 

For this sum, we found a slightly different weak change-of- variable formula for a process that acts similar 
to fBm with H — 1/4. In that case, the sum converges weakly to a stochastic integral plus a correction 
term in the form of an Ito integral of /", namely 

f(X t ) = f{X )+ [ f(X s )°dX s + l [ f"(X a )dB s , 
Jo z Jo 

where B is a scaled Brownian motion, independent of X, with a given variance. As is suggested in [5], 
there are other forms of Riemann sums that can be tried, including one for any fBm with H > 1/10. 
We expect that our theoretical tools could be applied to the H = 1/10 case as well, but this is not 
pursued in the present paper. 

A brief outline of the paper is as follows. In Section 2, we give some background on Malliavin 
calculus, and the main analytical tools that will be used in this paper. We also recall the central 
limit theorem (proved in [6]) that provides the main theoretical basis for our result. In Section 3, we 
identify the covariance conditions on the process X for the CLT to hold, and prove The proof 
essentially consists of restating the Riemann sum (jTJ) as a sequence of terms dominated in probability 
by 3-fold Skorohod integrals; then verifying that the CLT conditions are met. Section 4 discusses the 
three examples listed above, and demonstrates the procedure for verifying the covariance conditions. 
In Section 5 we give proofs for three of the longer lemmas from Section 3. 

The main inspirations for this paper were [5] and [TU]. Most of the notation follows that in [3] and 

®- 

2 Preliminaries and notation 

Let X = {X(t),t > 0} be a centered Gaussian process defined on a probability space (O, J 7 , P) with 
continuous covariance function 

E[X(t)X(s)\ = R(t,s), s,t > 0. 

We will always assume that T is the er— algebra generated by X. Let S denote the set of step functions 
on [0, T] for T > 0; and let f) be the Hilbert space defined as the closure of £ with respect to the scalar 
product 

(l[o,t],l[o, s ])« =R(t,s), s,t>0. 
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The mapping l[o,t] ► X(t) can be extended to a linear isometry between Sj and the Gaussian space 
spanned by W. We denote this isometry by h ^ X(h). In this way, {X(h),h £ ij} is an isonormal 
Gaussian process. For integers q > 1, let fj 09 denote the q th tensor product of Sj. We use f) 09 to denote 
the symmetric tensor product. 

For integers q> 1, let T-L q be the q th Wiener chaos of X, that is, the closed linear subspace of L 2 (Vt) 
generated by the random variables {H q (X(h)),h £ S),\\h\\fj = 1}, where H q (x) is the q th Hermite 
polynomial, defined as 

x 2 d q ^ 

In particular, H\{x) = x, Hi(x) — x 2 — 1, and H^(x) = a; 3 — 3a;. For q > 1, it is known that the map 

I q (h®«)=H q (X(h)) (4) 

provides an isometry between the symmetric product space .f) 09 (equipped with the modified norm 
^j=j | ■ \\s~,®q) and W q . By convention, Hq = R and Iq(x) = x. 

2.1 Elements of Malliavin Calculus 

Following is a brief description of some identities that will be used in the paper. The reader may refer 
to [9] for a brief survey, or to [TT] for detailed coverage of this topic. Let S be the set of all smooth 
and cylindrical random variables of the form F = g(X(<pi), . . . , X(4> n )), where n > 1; g : R™ — > E is an 
infinitely differentiable function with compact support, and <pi £ Sj. The Malliavin derivative of F with 
respect to X is the element of L 2 (il,f)) defined as 

n Pi 

DF = Y J -^{X{ ( j> 1 ),...,X{ct> n )) ( l )i . 

i—l 

In particular, DX(h) = h. By iteration, for any integer q > 1 we can define the q th derivative D q F, 
which is an element of L 2 (il, fj 09 ). 

For any integer q > 1 and real number p > 1, let 1D> 9 ' P denote the closure of S with respect to the 
norm || • ||b<j.p defined as 

||F||^=E[|^H+^E[||^F||^]. 

We denote by S the Skorohod integral, which is defined as the adjoint of the operator D. This operator 
is also referred to as the divergence operator in fTI]. A random element u £ L 2 (n,$j) belongs to the 
domain of <5, Dom <5, if and only if, 

\E[(DF,u)^] \ < c u ^E[F 2 ] 

for any F £ ED 1 ' 2 , where c u is a constant which depends only on u. If u £ Dom S, then the random 
variable S(u) £ L 2 (fl) is defined for all F £ D 1 ' 2 by the duality relationship, 

E[F5(u)] = E[{DF,u) Si ] . 

This is sometimes called the Malliavin integration by parts formula. We iteratively define the multiple 
Skorohod integral for q > 1 as S(S q ~ 1 (u)), with 6°(u) — u. For this definition we have, 

E[FS q (u)]=E[(D q Fu)^ q ], (5) 

where u £ Dom 8 q and F £ D 9 ' 2 . Moreover, if h £ f) 09 , then we have 8 q (h) = I q (h). 
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For / £ fj 0p and g £ Sj Qq , the following integral multiplication formula holds: 

pAq 



r=0 



s p (f)s q (g) = E H ( J) ( J )^ +9 " 2r (/ ®r 5), (6) 



where <S> r is the contraction operator (see, e.g., [llj . Sec. 1.1). 

We will use the Meyer inequality for the Skorohod integral, (see, for example Prop. 1.5.7 of |llj). 
Let D fc,p (io® fc ) denote the corresponding Sobolev space of .?)® fe -valued random variables. Then for p > 1 
and integers k > q> 1, we have, 

||<5 9 (m)||b*=-<3.p < Ck, P \\u\\ a k,p(jj® q ) (7) 

for all u £ D fe ' p (io lg ' fe ) and some constant Ck, P - 

The following three results are well known, and will be used extensively in this paper. The reader 
may refer to [9] and [11] for details. 

Lemma 2.1. Let q > 1 be an integer. 

(a) Assume F £ W' 2 , u is a symmetric element o/Dom 5 q , and (D r F, 5^ (u)) fli8r £ L 2 (Q,S)® q ~ r ~i) 
for all < r + j < q. Then (D r F, u) ^ 8r £ Dom 5 r and 

q 



r=0 ^ ' 

(b) Suppose that u is a symmetric element of W +k ' 2 (f)®^) . Then we have, 



i=0 



(c) Let u,v be symmetric functions in D 29 ' 2 ^® 9 ). Then 

q / \ 2 



i=0 

In particular, 



E[6 q (u)6 q (v)] = E (l) E [( D9 ~ iu > D9 ~ iv )z®(^ 



\S"(u)\\l Hn) =E[S"(u) 2 ]<Y^ 



1 / \ 2 

2 



i=0 



D q ~ l u\\ 



2.2 Criteria for convergence in the Skorohod space D[0,oo) 

Definition 2.2. Assume F n is a sequence of d— dimensional random variables defined on a probability 
space (f2, J 7 , P), and F is a d— dimensional random variable defined on (f2, P), where F C G- We say 
that F n converges stably to F as n — >■ oo, if, for any continuous and bounded function / : M d — > M and 
R-valued, J 7 — measurable random variable Z, we have 

lim E(f(F n )Z)=E(f(F)Z). 



The following central limit theorem first appeared in [5], and the present multi-dimensional version 
was proved in [5J. This will be the main theoretical tool of the paper. 
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Theorem 2.3. Let q > 1 be an integer, and suppose that F n is a sequence of random variables in R d 
of the form F n — 8 q {u n ) = (8 q (u n ), . . . , S q (u n )^j, for a sequence of M. d — valued symmetric functions u n 
in J} 2q ' 2q ($)® q ). Suppose that the sequence F n is bounded in L (QfS)) and that: 

(a) (vP n , &j£=i(D ae F-lf) ® hj ^ m converges to zero in i 1 (f2) for all integers 1 < j,j( < d, all integers 
1 < oi, . . . , a m , r < q — 1 such that a\ + • • • + a m + r = q; and all h £ $j® r . 

(b) For each 1 < i,j < d, (u l n ,D q F^j^ m converges in L 1 (f2,io) to a random variable sy, such that the 
matrix £ := (sij) dxd is nonnegative definite (that is, A T SA > for all nonzero A £ R d ). 

Then F n converges stably to a random variable in R d with conditional Gaussian law Af(0, S) given X. 

Remark 2.4. Conditions (a) and (b) mean that for q > 1, some combinations of lower-order derivative 
products are negligible. In particular, for q = 3, then the following scalar products will converge to zero 
in L 1 ^,^): 

. (K,h)^ 3 for all h £ 3® 3 . 

• (u n ,DF n ® for all h £ $j® 2 and all j (including i = j). 

• (u l n , D 2 F n ® h) ^ m for all /i € i5 and all j. 

. «,£> 2 i^ ® I?J*> a ., and (u n ,DF n ® #F„ fc ® Df*) a .,for all 1 < fc,j,*< d. 
Only the 3 rd -order derivative products (u l n , D 3 F-fy @ 3 converge to a nontrivial random variable. 

Remark 2.5. It suffices to impose condition (a) for h £ So, where So is a total subset of $j® r . 

For the main result of this paper, we will also require that the sequence F n satisfies a relative 
compactness condition in order to establish convergence in the Skorohod space D[0,oo). 

Corollary 2.6. Suppose {G n (t),t > 0} is a sequence of H-valued processes of the form G n (t) = 
5 q (u n (t)), where u n (t) is a sequence of symmetric functions in H> 2q ' 2q ($j® q ). Assume that for any finite 
set of times {Q = to < t\ < ■ ■ ■ < td}, the sequence 

(G n (ti) — G n (to), ■ ■ ■ , G n {td) — G n (td-i)) 

satisfies Theorem 2.3; where the dxd matrix X is diagonal with entries s (ti) — s (ti—\). Suppose further 
that there exist real numbers C > 0, 7 > 7 and ft > 1 such that for each n and for any < t\ < t < t%, 
we have 

E {\G n (t) - G n (h)r \G n (t 2 ) G n (t)f<] < C ( L " t2j = L " tlJ 

Then the family of stochastic processes {G n , n > 1} converges as n — > 00 to the process G = {Gt, t > 0}, 
where G(t) is a Gaussian random variable with mean zero and variance s (t). Equivalently, we can say 

that G n {t) -=+ V^t) Z as n — > 00, where Z ~ A/"(0, 1). 

This convergence criteria in D is well known (see, e.g, [2], Thm. 13.5). 
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3 Convergence of the Stratonovich integral 

3.1 Covariance conditions 

Consider a Gaussian stochastic process X :— {X t ,t > 0} with covariance function E [A S X(] = R(s,t). 
Assume R(s,t) satisfies the following bounds: for any T > 0, < s < 1, and s < r,t < T: 

(i) E [(X t — X t - S ) 2 ] < Cisi , for a positive constant C\. 

(ii) If t > s, 

\E[X^Xl s ]\<C 2S ^ e (t- S )- s 

for some C2 and 1/2 < 8 < 1. 

(iii) For £ > 4s, 

|E [(X t - X t _ s ) 2 - (X t _ s - AV 2s ) 2 ] I < C 3 s* +!/ (t - 2s)-" 
for some constants C3 and v > 1. 

(iv) There is a constant C4 and a real number A £ (i, i] such that 



\E[X r (X t -X t _ s )]\ < 



Ci-s ((< - s) A -! + |< - r| A - x ) if |t - r\ > 2s and t > 2s 
CiS X otherwise 



(v) There is a constant C5 and a real number 7 > 1 such that for t A r > 2s and \t — r\ > 2s, 

|E [(X t - X t _ s )(X r - X r _ s )]| < C 5 s* +7 |i - r|-T. 

(vi) For integers n > and integers < j, k < nT, define /3 n (j, k) := E (Xj+i_ — X j_)(X k+i_ — Xk_) 
Then for each real number < t < T, 

[ntj-l 

lim V p n (j,k) 3 =r,(t), (8) 
j,k=a 

where -q(t) is a continuous and nondecreasing function with 77(0) = 0. As we will see, r](t) is 
comparable to the 'cubic variation' [X,X,X] t discussed in [S] and [ID] . As described in [TU], these 
terms are related by Theorem 10 of |12) . 

In particular, it can be shown that the above conditions are satisfied by fBm with Hurst parameter 
H = 1/6. In Section 4 we will show additional examples. 

In addition to conditions (i) - (vi) on X, we will also assume the following condition (0) on the test 
function /: 

(0) Let / : K — )• M be a C°° function, such that / and all its derivatives have at most polynomial growth. 

Consider a uniform partition of [0,oo) with increment length 1/n. The Stratonovich integral of 
f'(W) will be defined as the limit in probability of the sequence (see |10j): 



LntJ-l 
3=0 



where AXj_ = Xj+i — Xj_ ■ 

The following is the major result of this section. 
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Theorem 3.1. Let f be a real function satisfying condition (0), and let X = {X tl t > 0} be a Gaussian 
process satisfying conditions (i) through (vi). Then: 

(X t , A (x t , f(X t ) - f(X ) + ^J* f^(Xs) dB^j 

as n — > oo in the Skorohod space D[Q,oo), where B = {Bt,t > 0} is a scaled Brownian motion, 
independent of X , and with variance E [Bf\ = rj(t) for the function r\ defined in condition (vi). 

The proof follows from Theorem 2.3 and Corollary 2.6, and is given in a series of lemmas. Following 
is an outline of the proof. After a preliminary technical lemma, we use a Taylor expansion to decompose 

$*(t) = f(X t ) - f(X ) + ±F n (t) + A„(t). 

We first show that A n (t) — > as n — > oo; then we show that F n (t) satisfies Theorem 2.3. Next we 
show that F n (t) + A„(t) is relatively compact in the sense of Corollary 2.6, and the result follows. 



Introduce the following notation, which is similar to that of [5] and [5]. Let e t — l[o,t]! di 



let Xj_ = i I, + Xj+i J ; and i± = \{s j_ + Sj+i ). In the following, the term C represents a generic 



positive constant, which may change from line to line. The constant C may depend on T and the 
constants in conditions (0) and (i) - (vi) listed above. By the isometry between the space generated by 
X and the Hilbcrt space fj, we will use the terms E[X s Xt] = (e s ,St)^ interchangeably 

We begin with the following technical results, which follow from conditions (i) through (v). 

Lemma 3.2. Assume {^,0 < t < T} satisfies conditions (i), (ii), (iv) and (v). For integers n > 1, 
r > 1 and integers 0<a<b<c< [nT\ , there exists a constant C > 0, which does not depend on 
a, b, c or r, such that: 



(a) 



(b) 



sup 

0<j,fc<[nTJ 



dj_,dk) 

r» / S 



< Gn 3 ; and sup sup 

0<u<T 0<j<[nT\ 



< Cn 



O 

<Cn-^{b-a + lf- 6 \ and 



b 



< Cn 3 for r > 1. 



(10) 
(11) 



3=a 



(c) ForO<u,v< T, 



(d) 



< C: and 



< Cr 



-2A 



3=a 



b 

\(dj_,d±) r < C(b-a+i)rri. 



j,k=a 



(12) 
(13) 

(14) 
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(e) 



(15) 



k=b+l j=a 



where e = max{l — 9,2 — 7}. 



Proof. We may assume a = 0. For part (a), the first inequality follows immediately from condition 
(i) and Cauchy-Schwarz; and the second inequality is just a restatement of condition (iv). For (b), 
applying condition (i) for j = and condition (ii) for j > 1, we have: 



EH 



j=o 



cb 



Then if r > 2, 



E | E [*2±i 



J=0 



^0 

< Cn-^b+l) 1 - 9 . 



<Cn-^^j- re +Cn 
3=1 



because 9 > 1/2 implies j~ r0 is summable. 

For (c), define the set J c — {j : < j < b, j — or \j — nu\ < 2 or \j — nv\ < 2}, and note that 
I J c | < 7. Then we have by (a) and condition (iv), 



E \{^ d i) I ^ E + Cn ~ x E ,A_1 + ii - 

< Cn- X + Cn- x (b+l) x < C, 



nu\ x - 1 ) 



and 



E ( e »' a *> \ £ "' 9 i) - E c^ 2A + C ^ 2A E 0' 2A ~ 2 + Ij - h 2A ~ 2 + b' - h 2A ~ 2 ) 

00 

P =i 

< Cn- 2A 

because A < 1/3. 

For (d), define the set: Jd — {j, k : j A k < 1 or |j — fc| < 2} , and note that | J d \ < 6(b + 1). Then 
we have by (a) and condition (v) 



E \\ d i ,d % 

j,k=0 



< sup 

3-k 



dj_,dk 



3,k=0 



j,k=0 



<Cn-^ ^ \j - k\~~< 

\ti,k)£J d (j,k)<fJ d 

< C(6 + l)n~S. 
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In particular, if r = 3 and b = \nt\ — 1 (as in condition (vi)), the sum converges absolutely, and the 
sum vanishes if r > 3. 

For (e), we consider the maximal case, which occurs when a = 0: 



b-l 



e ek^ 

k=b+l j=0 



E E + E 



E 

fe=&+l k=b+l j=l ' "' k=b+l 

Note that do = e±. By part (c) and condition (v), respectively, this is 

c 6-1 



< E 



k=b+l k=b+l j = l k=b+l 

< Cn-^{c - b) 1 - + Cn~5 ( c _ 6) 2 -7 + c n ~i 

< Cn-i{c-b) e , 

where e — max{l — 9, 2 — 7} < 1. □ 
3.2 Taylor expansion of 

The details of this expansion were mainly inspired by Lemma 5.2 of [10] . We begin with the telescoping 
series, 

f(x t ) = /(o) + f(x t ) - + [f( x ^ - f( x i) ■ 

By continuity of / and X, we know that for large n, f(X t ) — f(X y n t\ ) — > uniformly on compacts in 
probability (ucp), so this term may be neglected. For each j, we use a Taylor expansion of order 6 with 



residual term. Let hj 



X j+i — X j_ 



Then: 



/(*i±i) - f{Xi) = [f{Xi. + hj) - /(Jfi)J - - h 3 ) - f(Xj_ 



h : 



E/ (fc) (^)^ + ^01 - ( E(- 1 ) fc / (fc) (^)-^ + ^t(j) 



k=l 



\k=l 



= fiX^AXj, + l/^X^AXj + J^/W(^)AXi + - 

where AJi = 2/ij; and R^(j), R^(j) are Taylor series remainder terms of order 7. Next we compute: 

f'(X i±1 ) + f{X L ) . 1 / - . \ 1 / - - \ 
=-2 =- - f(Xi) = - (/'(Xj. + hi) + - (/'(Xi - ft,) - f(Xi)) 

= |E/ (1+fc) (^)f +^(i) + ^E(-i)V (1+fc) (^)| +K~U) 



k=l 



k=l 



= i/( 3 )(X,)AX| + J^/(5) ( X 4 )AX| + i (A+ + KJ) , 
where Kf,K~ are remainder terms of order 6. Combining the two equations, we obtain 



/(X i±i )-/(X i ) 



/'(Xi±i) + f'{X ± ) 



^AX^ - l/( 3 )(X,)AXi - JL/(5)(^)AXi 



12 

- «n(i) - \ [k+U) + k-(j)] AXx. 



Our first task is to show that the terms and the remainder term vanish in probability. 
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Lemma 3.3. For each integer n>\ and real numbers < ti < t 2 < T, 

2" 



E 



\nt 2 \-\ 

J2 / (5) (^)AX| 
vj=L nt iJ 



< Cn~3 ([nt 2 \ - [ntij). 



(16) 



The proof of this lemma is technical, and is deferred to Section 5. 



Lemma 3.4. For integers n > 1, let 



[nt]-l 

Zn(t)= 

3=0 



Ri(j) - + 2 (K+ti) + K~{j)) AX, 



Then for real numbers < t x < t 2 < T, we have 



E 



(17) 



(Z„(i 2 ) - Z„(ti)rj < Cn—z (\nt 2 \ - \nt Y \Y . 

Proof. We may assume ti = 0. Observe that each term in the sum Z n (t) has the form 

Cf^^)AX\, 

where £j is an intermediate value between X± and Xj+i . Using the Holder inequality, for each < 
j, k < [nt 2 \ we have 



E 



/ (7) te)/ (7) (&)AXlAX! 



< ( sup E f {1) {uX L + (1 - ?i)X 2 +i) 4 sup E f [7) (vX K + (1 - v)X h±1 f 

V0<M<1 In J Q<-U<1 L 



E 



AXl 



E 



AXl 



By condition (0), the first two terms are bounded. By condition (i), E 
by the Gaussian moments formula that 



AXl 



<C\n 3 ; and we have 



E 



AX 



2<S 



/ i\14 

<27!!(Cin-3j , 



hence it follows that 



[nt 2 \-l [nt 2 \-1 

C £ E[/( 7 )fe)/ (7) (6)AXlAX|] <C n-i <C[nt 2 \ 2 n-i. 

j,k=0 j,k=0 



□ 



3.3 Malliavin calculus representation of 3 rd order term 

From Lemmas 3.3 and 3.4, it follows that as n becomes large, then f(X t ) behaves asymptotically as 



\nt\-l 



\nt\-l 



3=0 



3=0 



We now turn to the f^(Xj_) term in the Taylor expansion. 
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Remark 3.5. It may happen that the upper bound of condition (v) is such that 



[nt]-l 

r,(t) < lim V |/?„(j, fc) 3 | =0 

n— too * — ' 1 1 

j,k=0 



for all i, which implies 



lim E 



|nij-l 

£ /( 3 )(^)ax| 



for any function / satisfying condition (0). In this case, converges in probability and we have the 

change-of- variable formula ([2]) . Indeed, this corresponds to the case of zero cubic variation discussed in 
[5]. In the rest of this section, we will assume that r/(t) is non-trivial. 

Consider the 3 rd Hermite polynomial H^{y) = y 3 — 3y. For y — \ \ax\ \ 2 ' ^ follows that 

\nt\-l \nt\-X / /\X j 



ratj-l L«*J— 1 / 

£ /( 3 >(A^)AAT 3 = £ ||AXi||| 2 /( 3 )(X 4 )if 3 i 

7=0 7=0 V 



|nt|-l 

+ 3 £ \\&Xj_f L2 fW(Xj_)AX i 

3=0 

The second term is dealt with in the next lemma. The proof is technical, and is deferred to Section 5. 
Lemma 3.6. For integers n > 1 and integers < a < b < nT, 

6-1 



E 



£||AX i ||| 2 /( 3 )(X i )AX J 



< Cn-s. 



Next, we consider the term. By (0| and Lemma 2.1. a we have 



LntJ-l / AX 

]T WAX^hf^iXM' 

3=0 



\AXo 



3=0 

LntJ-1 - 

£ <S 3 (/ (3) (Xi)3f) +35 2 ((£>/( 3 )(Xi),9f) 
+ 3^( J D 2 / (3) (Xi),9f )^ 2 S j+(D 3 fW(X i ),df ) 

ntj-l 

E * 3 (/ (3) (^i)flf 3 )+^n(*). 



3=0 



As n — >■ oo, we show that the term P n (t) vanishes in probability. 
Lemma 3.7. For integers n > 1 and reaZ numbers § < t\ < t2 < T , 

E[P„(t) 2 ] <C(Lrrf 2 J - ln*iJ)n"*- 



Proof. We may assume t\ = 0. We want to show 



E 



E 



|nt 2 J-l 
3=0 



J=0 



and 



E 



|nt 2 J-l 

E (^ 3 / (3) (^),5f 



< C\nt 2 In- 



Proof of (|T8]). By p) we have 



E 



E (^/ (3) (^)' 



< E 



2E 



3®3 



L"*2j-1 , 

E ((^/ (3) (^),5f) ,(l?/C 3 >(%,9 

|nt 2 J-l , 

E ( (^ 2 / (3) (*i ) , flf 3 ) _ , (^/ (3) (*a ) , 5f 3 



_<yS2 



■E 



|nt 2 J-l 

E (£> 3 / (3) (^i),5f 

j,fc=0 



D 3 f( 3 \x h ),df 3 



< sup E 

o<j<[nt a } 



sup E 

Q<j<[nt 2 ] 



sup E 

0<j"<Lnt 2 J 



/ (4) (^) 2 " 



|nt 2 J-l 



sup {£j_,dj_) ) \d±,dk 

0<j<[nt 2 l \ » »/« ^ \ « » 



sup (ij_,d^) E (^i)^i 

0<j<|n* 2 J X " 1 * ' 



j,fc=0 



L"* 2 J-i 

/ (6) (% 2 ] E 



n I ir, 



By condition (0) and Lemma 3. 2. a and 3.2.d, 



sup E 

0<j<[nt 2 } 



sup E 

o<j<L"* 2 J 



2 L™*2j-1 2 

sup {£±,d±\ E \9±,dk) < Cn~^ \nt%\ ; and 
o<j<Lnt 2 J x ™ n/ a , „ \ n n 'Sj 



/ (5) (t) 2 sup (fi,5 4 > E 



4 LnfeJ-l 

a 



< Cn~3 [ n t 2 \. 
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Then by condition (0), Lemma 3. 2. a and 3.2.b, 



sup E 

0<j<jntaj 



\nt 2 \-l 

/ (6) (^) 2 ] E 

j,k=0 



i±,d± 

n n/ Sj 



£ k , d k 



< C sup (£±,dj_ 

0<j<\ntz\ x " ™ 1 * 



Ln*2j-1 
j=0 



Proof of ([TO)) and ([2H)l . The same estimates apply for the other terms, since 

2" 

/ LntaJ-1 \ \ 



E 



\nt 2 \-l 

E (^ 2 / (3) (ix),9f 

3=0 



< 



E 



L"*2j-1 

E 

j,fe=0 
Lni 2 J-l 



7 (3) (^),sf) ,(^ 2 / (3) (^),5? 3 ) 



E 



E (£> 3 / (3) (^),9f) 88 -(0 a / (3) (*i),0f 3 ) 



and (|20l) is bounded in the above computation as well. □ 

3.4 Weak convergence of non-trivial part of 3 rd order term 

We are now ready to apply Theorem 2.3 to the term 

[ntj— l 

E ^ 3 (/ (3) (^)9f 

3=0 

Let = to < t\ < ■ ■ ■ < td < T be a finite set of real numbers, and for i = 1, . . . , d define 

[ntij-i 

<= E / (3) (^)9f; 

j=[nti-ij 

and define the d— dimensional vector F n — (F^, . . . , F^), where each — S 3 (u l n ). 

To satisfy the CLT conditions, we must deal with terms of the following forms (see Remark 2.4): 

1. (K,h)^ 3 iorhe^ 3 , 

2. {ui,DFi®h)^ 3 for/te^ 2 , 

3. (ui,D 2 Fj <g> hi) a9a + (K,DFi ® #F„ fc <g> & 2 ) fl „, for h u h 2 e 0, and 
4 <V D 3 F' 1 ') + <V D 3 F :i ') 

+ (<, # 2 ^ ® £F„ fc )^ 3 + «, ® DF„ fc ® DFi)^ 3 . 

We must show that all terms converge to zero except for the terms {u l n , £> 3 F^) J . 03 , i = 1, . . . , d, which 
will converge stably to a Gaussian random vector (Lemma 3.11). 
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Lemma 3.8. For each i,j,k,£= 1, . . ., d, the following terms vanish in i 1 (fl) as n — > oo; 

(a) (<, h)^ 3 for each h € f)® 3 . 

(b) (<, DFl ® h)^ 3 for each he?)® 2 . 

/ or h e Sj. 

Proof. We begin with two estimates that will be needed. For each 1 < i < d, 

E||.OF*||* < C; and 
E\\D 2 F^\\l m <C. 
Proof of (|2"Tj) . Let = [^ti-ij and 6i = [nti\ . By Lemma 2.1.b, 

DF; = <5 a (I)<)+3<5 2 (<). 

Hence, using Lemma 2.1.C, 

bi-l 



(21) 
(22) 



E 



|i^|| 2 <2 J] E |^<5 3 |/' 4 '(ii)3fj <5 3 

+ 18 E E [<5 2 (/( 3 )(li)9f )^ 2 (/ (3) (^)5f 



dj_,dk 



E 



bi-l 3 

j,k=a,i l=Q 

bi-l 2 

+ EE 

< c 

by condition (0) and Lemma 3.2.d. The proof of (|22|) follows the same lines, using Lemma 2.1.b to 
obtain 

D 2 F l n = <5 3 (£> 2 <) + 65 2 (£><) + 65«). 

Now for the main proof. Without loss of generality, we may assume that each h € Sj is of the form 
e T for some < r < T (see Remark 2.5). Then for (a) we have: 



E 



(<,h) 



E 



bi-l 



E (f {3) (X^)df,e T ®s u ®e V/ 

n—ai 

< sup {E f (3 \X^) } sup J2 |<^%), 

ai<m<bi ^ J T.m ^ J „,_„ ' 



bi-l 



-2X 



< Cn 

where we used Lemma 3. 2. a and Lemma 3.2.c. For (b), 



(<,DFi®Er®s a ) s .® 3 



< J sup E 

y ai<m<bi 



bi-l 



/< 3) (%) J2 (^(dT,DF^®e T ®e u 



m—ai 
bi-l 



<Csup||%|| sJE\\DFi\\% E <%>^M%' e «>i 



15 



By condition (i), \\dm \\fi < Cn e , and so by (f2"Tj) and Lemma 3.2.C we have an upper bound of Cn 
For (c), by similar reasoning along with Lemma 3.2.C and (|22|) . 



-i-2X 



bi-l 



E 



< J sup E 

y ai<m<bi 



6<-l 



<Csup||fe|| JE\\D*FM\% m T \(d^,e r ) 

m I 

rn—ai 

< Cn-*. 

The estimate is similar for the term E (u^, DF^ <E> DF% ® e t)^,{ 
Now we focus on the terms (u^, D 3 F£)^, S3 . By Lemma 2.1.b, 

D 3 FZ = S 3 (D 3 ui) + 95 2 (D 2 ui) + m(Du{) + Qui 
so that (u % n , D 3 F£)^ 3 can be written as, 



□ 



l.m 



+ 9 / (3) (*» (/ (s) )df) , i t ) (d^ , si ) 

+ 18^/( 3 )(^)<5(/( 4 )(^)5a) (%,%) (%^%) 2 +6«,<) 



:= An(i,j) + B„(«,i) + C n (i,i) +6(<,<) fli83 , 
where are the indices for u l n ,D 3 Fl, respectively, with [fiii-ij < m < \nti\, and [fiij-ij < I < 



Lemma 3.9. For each 1 < i, j < d we have 







as n — > oo, where 5ij is the Kronecker delta. 
Proof. We will show that for each 1 < i,j < d 

Urn E \A n (i,j)\ = lim E \B n (i,j)\ = lim E \C n (i,j)\ = 0. 



and moreover, if i ^ j then lim„ 



E 



("in 



0. 



To begin with, observe that if 5(2;) is a function satisfying condition (0), then it follows from condition 
(i) and Lemma 2.1.C that that for q = 1, 2, 3, 



sup 



6* (fl(*i)0f) 



< Csup 

L 2 J 



< Cn-i. 



(23) 



For the terms A n (i,j), B n (i,j), C n (i,j) we include the case i = j. We have 



E\A n (i,j)\< S up / (3) (lm) 



L- 



sup 



S 3 (f^iX^df) 



L 2 



sup 
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Using condition (0), (j2"3")l . and Lemma 3. 2. a, respectively, we have 



2 

) 

a 

so that Lemma 3.2.C gives E\A n (i,j)\ < Cn~2~ 3X . Next, using condition (0), (|23|) and Lemma 3. 2. a, 

2 

a ■ 



E\B n {i,j)\ <9sup /( 3 )(1 



sup 



(/( 5 )(Xi)9| 



sup < 

L 2 m \ 



<Cn-3- 2A ^|(^,a. 



and so by Lemma 3.2.d, 



E\B n (i,j)\ < Cn~^~^max{[nti\, \ntj\} 



which converges to zero since 2A > 1/3. Similarly for C n (i,j) using Lemma 3.2.d, 



E\C n (i,j)\ < 18 sup / (3) (^) 



L- 



sup 



sup 

L 2 Lm 



< Cn-i- x ^2^di.,d^^ < Crri~ x max{ [nU\ , [ntjj } < Cn 



E 

£,m 



For the second part, we may assume i < j. Using Lemma 3.2.e 

2 
L 



E 



<<»<) 



< sup 



E 



which converges to zero because e < 1. 

Lemma 3.10. Using the above notation, for each 1 < i,j,k,l < d we have 

l2 7711 o 7-1 77.fc\ 2 



lim E 



lim 



(ui,DFi®DF*®DF*) 



= 0, and 
0. 



2 

ij®3 



The proof of this lemma is deferred to Section 5. 



□ 

(24) 
(25) 



Lemmas 3.8, 3.9 and 3.10 show that condition (a) of Theorem 2.3 is satisfied, and moreover that 
the only non-trivial terms are of the form Q (u l n ,u l n ^ It remains to establish the convergence of 
these terms to a non- negative random variable 6s|. With this result, it follows from Theorem 2.3 that 
the couple (X, F n ) converges stably to (X, (), where C = (Ci) • • • ) Crf) is a vector whose components are 
conditionally independent Gaussian random variables with mean zero and variance 6s 2 . 

Lemma 3.11. For each 1 < i < d, conditioned on X , 



where each sf has the form 



lim (<,<) 



s(t 4 ) 2 - s(fi-i) 2 



f^(Xs)Mds). 
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It follows that on the subinterval [it— i,t»] we have the conditional result 

almost surely as n — > oo, which implies 

Fi-^V6 f % f {3) {X s )dB s , (26) 
•/ti-i 

where {B tl t >0} is a Brownian motion, independent of X, with variance 77(f). 



Proof. Let a = [rati- 1 J and 6 = [nt^ , and recall /3 n (J, k) = (^dj_,dkj , from condition (vi). We have 

6-1 

«,<)^ 3 = E / (3) (*i)/ (a) (*i)/3na*) a - 

j,fc=a 

For each n, define a discrete measure on {1, 2, . . . }® 2 by 

00 

3,*=0 



where denotes the Kronecker delta. It follows from condition (vi) that for each f > 0, 

/j([0,f] 2 ) := lim Mnflnt], LntJ) 



n— J- 00 

LntJ-l 



= lim V /3„(i,fc) 3 = 77(f). 

n * — ' 

j,fc=0 



Moreover, if < s < f then 



[nsj-l L«*J-i 

/j„ ([nsj, |nf|) = \i n ([ns\, [ns\) + ^ E Mj,k) 3 , 

j=0 k=[ns\ 

which converges to zero because the disjoint sum vanishes by Lemma 3.2.e. Hence we can conclude 
that fi n converges weakly to the measure given by /i([0, s) x [0, f]) = n(s At). It follows by continuity of 
f( 3 '(X t ) and Portmanteau theorem that 

IMI-1 , 
J2 / (3) (Xi)/ (3) (X,)^ n (j,fc) 3 = / f^(X s )f^(X u )l s<t l u<t fi n (ds,du) 
j,k=o " " • /r2 

converges in L 1 (i7,i5) to 



It follows that on the subinterval [i*_ i,ti] we have the result 

«,o^ 3 -^ r f {3) (x S ) 2 V (d S ) 

in L 1 (fi,^) as n — > 00. Using the Ito isometry for the above integral, we conclude (|2"6")l . □ 
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3.5 Relative compactness of the sequence F n (t) + A„(t) 
Let 

1 L " tJ_1 / x 



3=0 



To establish convergence of G n (t) in D[0, oo), we need to show that {G n (t)} is relatively compact. For 
this, it is enough to show that there exist real numbers a > 0, j3 > 1 such that for each T > and any 
< ti < t < t 2 < T we have, 







E [\G n (t) - G n (h)\ a |G„(t 2 ) - G n (t)\ a ] < C ( [nt2j - [nhl 
We will do this in several parts. From the preceding section we have, 

LntJ-l |n*J-l 

G„(t) = /(X ifI )-/(X ) + - ( 5 3 (/( 3 )(X,)af)+— £ f^iX^AX 

j=o " ' j=0 

Q L"-*J — 1 2 



J=0 



L 2 ' 

[ntj-l 



= /(x i ^)-/(x ) + l E <* 3 (/ (3) (* 



3=0 



A n (t). 



By Lemmas 3.3, 3.4, and 3.7 we have 



E 



'[nt 2 \-l 

E / (5) (^)AX| 
u'=L«*iJ 



< Cn-i ([nt 2 \ - [nti\); 



E [(Z n (t 2 ) - ^(tx)) 2 ] < Cn~l {[nt 2 \ - [nt^f ■ and 
E[P„(i) 2 ] <L7n-t([ni 2 J - [ntij). 

Each of these estimates has the form 



E 



(U n (t 2 ) UM) 2 } < On-? (|nt 2 J - L^iJ) f < C ( ^ KlJ 



where £ < /3 and /3 > 1, hence it follows by Cauchy-Schwarz that for £i < t < t 2 we have 



E 



[\U n (t) - U n (h)\ \U n (t 2 ) - U n (t)\] < C ( [nt2j - L "* lj y 



so each of these individual sequences is relatively compact. For the term, 



[nt\ - 1 

r„(t)= E \ Ax i 

3=0 



f {i) {X ± )5{d L ), 



L- 



we have by Lemma 3.6 that Y n (t) vanishes in probability. However, to show relative compactness we 
need a different estimate. 
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Lemma 3.12. For < t\ < t 2 < T such that \nt 2 \ — \ nti\ > I, we have 



E 



(r„(t 2 )-r„(t 1 )) 4 ] <Cn- 2 ([nt 2 \ - [nhlf + Cn-i-^dnhl - LniiJ)3 +4A . 
It follows that the sequence {Y n (t)} is relatively compact. 

Proof. Let $„ := $„(ji, j 2 , J3, ji) = I\t=i / (3) (^Zi), and let a = L n *iJ> b = L^J- We have 

4" 



E 



' b-1 



J2\\ AX i L2 / (3) (^w%) 



< sup AX 2 

a<j<b" n|lL2 



^ |e [Mh,h,h,j4)5(d^)5(d^)S(d^)S(d^) 



31,32,33,33 



< 



< 



Cn-s ^ |E[(£>[$ n 5(0£)5(0£)J(0£) 

jl ,32,33,3* 

4 

S H| E [ $ « ),5 (%) ,5 (%) (5 (%) 

jl ,32,33,34 r=l 

+ 3Cn"i ^ |e [$ n <K%)<K%) 



31,32,33,34 
4 



E E| E [( jD K r),5( % )<5( % ) 



31,32,33,34 r=l 



+ 3Cn-^ ^ |e[(d[$„5(0£) 



dj3,du 



£ jr. , 9i4 



where 



^ r) =/ (4) (^)f[/ (3) (^)- 



Continuing this process, we obtain terms of the form: 



31,32,33 ,34 



n I Sj \ n n / ft 

Cn~i ^ |E[a 2 $] (ii±,dii) U^dto) (d^di 



, and 



31,32,33,34 



^ | E [a 4 $] (i^,dA U^di,) U^d^) U^d 



31,32,33,34 



where <9 fc <I> represents the appropriate fc*' 1 derivative of By Lemma 3.2.C and 3.2.d, the sums of each 
type have, respectively, upper bounds of the form 

Cn- 2 (b - a)- 2 + Cn-i{b - a) + Cn-^ 4A (6 - a) 4X , 

hence we conclude that 



E 



(r„(t 2 )-y„(i!)) 4 ] < C{n- 2 ([nt 2 \ - LniiJ) 2 +n-t- 4A (Lnt 2 J - NiJ) l+4A ) , 



As for above terms, it follows by Cauchy-Schwarz that 



E 



Y n (t) - Y n {h)\ 2 \Y n (t 2 ) - Y n (t)\ 2 ] < C (n- 2 ([nt 2 \ - [nt^) 2 + n --^ AX {\_nt 2 \ - Ln^J)i+ 4A ) 



and thus {Y n (t)} is relatively compact. 



□ 
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Tightness of F n . 

To conclude the proof of Theorem 3.1, we want to show that the sequence {F n (t)} satisfies the 
relative compactness condition. 

Lemma 3.13. For < t x < t 2 < T, write 

F n (t 2 ) - F n (h) = J2 53 {f {3) (*^ d T) 

j=|ntij+l 

Then given 0<ti<t<t 2 <T, there exists a positive constant C such that 

E [\F n (t) - f»fti)| 2 |F»fe) - F n (t)\ 2 ] < C ( [nh{ = [n * lJ V . (27) 

Proof. We begin with a general claim about the norm of DF n . Suppose a, b are nonnegative integers. 
Let 

Lnt 2 J-l 

g a = J2 f {a) (X±)df. 
i=L«*iJ 

Then we have 

n\D b 9a\\%, 3+ A<c[ [n ^l [n ^ )\ (28) 
Proof of (f2"5)) . For each b we can write 



E 



\D b g a \\'f ) ®3+b) 



|n* 2 J-l 

£ 



< sup 

L«tiJ<j<L"*2j 



E 



4\ 2 



sup 



£ i. , £_fc 



26 



L«t 2 J-i 



< Cn- 2 {\nt 2 \ - L«*lJ)' 



by Lemma 3.2.d. 

Proof of ([27|) . By the Meyer inequality ([7]) there exists a constant C24 such that 



E 



where in this case, 



(<5 3 (W„)) 4 < C 3: 4||u„||„3,4( fl ®3), 

Un = £ /^(x,)af 

j=L«*iJ 



and 



E\\D z u 
,2 



2„, 1 1 4 



+ E\\Du n \\ 



l-D"' 



From d28j) we have E||£» b M„||4 03+i) < C h n~ 2 {\nt 2 \ - |n*ij)* for b = 0,1,2,3. From this result, given 
0<<i < t < t 2 , it follows from the Holder inequality that 



E \\F n (t) - F n {ti)\ 2 \F n {t 2 ) - F n (t)\ 2 } < C 



[nt 2 \ - [nti\ 



□ 
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4 Examples of suitable processes 
4.1 Bifractional Brownian motion 

The bifractional Brownian motion is a generalization of fractional Brownian motion, first introduced 
by Houdre and Villa [7J. It is defined as a centered Gaussian process B H ' K = {B?' K ,t > 0},with 
covariance given by, 

where H £ (0, 1), K e (0, 1] (Note that the case K = 1 corresponds to fractional Brownian motion with 
Hurst parameter H). The reader may refer to [13] and [8] for further discussion of properties. 

In this section, we show that the results of Section 3 are valid for bifractional Brownian motion with 
parameter values H, K such that HK > 1/6. 

Proposition 4.1. Let B t = {E>f' K ,t > 0} be a bifractional Brownian motion with parameters H,K 
satisfying HK = 1/6. Then conditions (i) - (v) are satisfied, with 9 — 2/3; A = 1/3; 

[5/3 if H< 1/2 , [2/Z + 2H if H < 1/2 and K < 1 

v = < , ,; and 7 = < 

if H> 1/2 [5/3 otherwise 

Proof. Condition (i). From Prop. 3.1 of [7J we have 

E [(£? f - B t - S ) 2 ] <Cs 2HK =Csi. 

Condition (ii). By Fundamental Theorem of Calculus, 

|E [B 2 - B 2 _ s ] I = t 2HK - (t - S ) 2HK = T 2^X(f + C) 2HK - X ^ < Cs(t - 

•/ — s 

Condition (Hi). 

E - 5 t _ s ) 2 - (£ t _ s - B t _ 2s ) 2 ] = E [{B t - B t . 2s )(B t - 2B t _ s + B t _ 2s )} 
= t 2HK i V 2 " + (* " *H K + i + (* " 2s) 2 "] * 



i [t™ + (< - 2sf H ] K + ± [(* - ,s) 2 * + (t - 2 S ) 2 «] * . 



In absolute value, this is bounded by 



1 



[t 2H + t 2H ] K - 2 [t 2H + (t - s) 2ff ] * + [t 2H + (t- 2s) 2H ] K 



1 



[t 2H + (t- 2s) 2H ] K -2[(t- s) 2H + (t- 2s) 2H ] K + [(t - 2s) 2H + (t- 2s) 2H ] K 



2 h 

Both terms have the form 

2- K \g(t)-2g(t-s)+g(t-2s)\ <Cs 2 sup \g"(x)\, 

xe[t-2s,t] 

for g(x) = (c + x 2H ) K . We show an upper bound for the first term g"{x), with the other one similar. 
We have 

sup \g"{x)\ < 4H 2 K(1 - K) [t 2H + x 2H ] K ~ 2 x 4 "' 2 + 2HK \2H - 1\ [t 2H + x 2H ] K ~ l x 2 "' 2 . 

x£[t-2s,t] 
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For the above values of H, K we have 



sup 2HK \2H-1\ [t 2H + x 2iI Y Y x 2H - 2 < C(t - 2s) 2HK - 2 . 

x£[t-2s,t] 

For the first term, if H < 1/2 then 

sup 4H 2 K(1 - K) [t 2H + x 2H ] K ~ 2 x 4H - 2 < C(t - 2s ) 2HK - 4H+4H - 2 
xe[t-2s,t] 

= C(t-2s) 2HK - 2 . 
On the other hand, if H > 1/2, then t > 4s implies t > 2(t — 2s), hence 

sup 4H 2 K(1 - K) [t 2H + x 2H ] K ~ 2 x 4H - 2 < 4H 2 K(1 - K)3 K - 2 (t - 2s) 2H ^ K - 2 K 4H - 2 



x 



xe[t-2s,t] 

< C(t - 2s) 2HK - 4H = C(t - 2s)^ 4H . 

Condition (iv). First, for the case |t — r"| < 2s or t < 2s, we have 
\E[B r (B t -B t _ s )}\ = \E [B r B t - B r B t - s ]\ 

< ± ( [r 2H + t 2H ] K [r 2H + (t- s) 2H ] K )+^k \V-t + s\ 2HK - \r-t\ 2HK \ 

< Cs 2HK - Cs? 

using the inequality a r — b r < (a — b) r for < r < 1. For \t — r\ > 2s, t > 2s, we consider two cases. 
First, assume r > t + 2s. 



\ 2HK -\r-t\ 2HK \ 



|E [B r (B t - B t -,)]\ = |E [B r B t - B r B t - s }\ 

< ^ ( V 2H + n K V 2H + (* - ')"*] *) + ^ 1 1'" * + < 
= ± 2HK [r 2H + (t + 2H ] ^ (t + t) 2 *- 1 ** + -L jf 2HK(r t + V f HK ~ 

< 2 1 - K HKs{t - s)-i + 2 1 - K HKs(r - 

where we used the fact that r — t > 2s implies r — t > 2(r — t — s). On the other hand, if r < t — 2s, 
then the estimate for 

^([r 2H +t 2H ] K -[r 2H + (t-s) 2H ] K ) 



is the same, and for the other term we have, 
1 



i I \r - t + s\ 2HK -\ r -t\* HK \<±f J° 2HK(t - r - ^ HK ~ 1 d^ 



< 2 1 - K HKs(t - r - s) 2 ™- 1 < 2i~ K HKs{t - r)"*, 

hence for either case we have an upper bound of Cs (if — s) A_1 + \t — r| A_1 ) for A = i. 
Condition (v). Assume t Ar > 2s and \t — r\ > 2s. We have 

E[(B t - B t _ s )(B r - B r _ s )\ 

-L ( [t 2H + r 2H } K [t 2H + (r- s) 2H ] K [{t s) 2H + r 2H } K + [(t - s) 2H + (r - s) 2H ] «) 



2 K 

1 

2 K 



+ 4 (\t- r + s\ 2HK - 2\t - r\ 2HK + \t-r- s\ 2HK ) . 
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This can be interpreted as the sum of a position term, ^-(p(t, r, s), and a distance term, ^-ipit — r, s), 
where 

<p(t, r, s) = [t 2H + r 2H ] K - [t 2H + (r - s) 2H ] K - [(t - s) 2H + r 2H ] K + [(t - s) 2H + (r - s) 2H ] * ; and 

1>(t -r,s) = \t-r+ s\ 2HK - 2\t - r\ 2HK + \t-r- s\ 2HK . 

We begin with the position term. Note that if K = 1, then ip(t, r, s) — 0, so we may assume K < 1 and 
H > g. Without loss of generality, assume < 2s < r < t. We can write </?(£, r, s) as 



2#if 



o 



([i 2 * + (r- 2H ] K_1 (r - fl 2 "" 1 - [(« - s? H + (r- if 11 ]*- 1 (r fl 3 *" 1 ) « 
4ff 2 *T(l - X) [(f - t?) 2 " + (r - f )2^] K - 2 { t - v ) 2H -\ r - ^ 2H ~ 1 <% dr,, 

o Jo 
so that 

\<p(t,r,s)\ < m 2 K{\ - K)s 2 [(t - s) 2H + (r- s) 2H ] K ~ 2 (t - s) 2H -\r - a) 2 *" 1 . (29) 
Using (|29l) . there are 3 cases to consider: 

• If H < 1/2, then for 2s < r < t — 2s, we have t — r < t — s and 

CS 2 [(i - s) 2H + (r - S) 2H ] K " 2 (t - S) 2 "- 1 ^ - S) 2 ^ 1 < Cs 2 (t - rfHK-IH-l^ _ s) 2H-x 

= C( — ) s ^+i +2ff (t-r)-i- 2H 



r — s 



■ r 



where 7 = | + 2H > 1. 



3 

• If iT = 1/2, then X = 1/3 and for 2s < r < t - 2s 

s 2 [(t - s) 2H + (r - s) 2ff ]^ 2 (i - S ) 2H -\r - sf H ~ l < s 2 \t - r|"*. 

• If H > 1/2, then note that for 2s < r < t - 2s 

s 2 [(t - s) 2H + (r - s) 2H ] K ~ 2 (t - s) 2H -\r - s) 2 ^ 1 < s 2 (t - s) 2HK - 2 < s 2 \t - r|~*. 

Next, consider the distance term ip(t — r, s). Without loss of generality, assume 2s < r < t — 2s. We 
have 

m - r ,s)\ = \\t-r+ s\ 2HK - 2\t - r\ 2HK + \t-r- s\ 2HK \ 

I ^ 2HK{2HK -l)[t-r + ti] 2HK ~ 2 dr\ d£ 

< Cs 2 (t - r - s) 2HK - 2 < Cs 2 \t - rj-i, 

since \t — r\ > 2s implies (t — r — s) _ § < 2§\t — r|~t. Note that when K < 1, then H < 1/2 implies 
7 < 5/3, so the upper bound is controlled by cp(t, r, s) in this K = 1 case. □ 
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Proposition 4.2. Let {B^' K ,t > 0} be a bifractional Brownian motion with parameters if < 1/2 and 
HK = 1/6. TTien Condition (vi) of Section 4 holds, with the function rj{t) = Cxt, where 



C K = 8 + 2 (( m + - +{m- 1)* 

\ m— 1 



Proof. First of all, we write 



[ntj-l J — l L rl *J — 1 

£ /3„(j,fc) 3 = 2 £ /3„(j,0) 3 + £ /3„(j,fc) 3 

j,fe=0 j=0 j,fe=l 



When j > 2, we have 



< 

" 2^ 



^ ([i + (i + iH A " - [i + | (?■ - D 2 "* - f Hh 

l — [\ H K[l + (j + xf H ] K - X (j + xf H -Hx +-^f 2HK(j 1 + yf HK -Hy 



< Cn~^{j - 1)"T. 
Therefore, using Lemma 3. 2. a for /3„(0,0) and /3 n (l,0), 



[ntj-l L nt J _1 

\Pn(j,0) 3 \ <2Cn- 1 + Cn- 1 ^ - I)" 2 < Cn" 1 ; 

3=0 i=2 



and in the rest of the proof we will always assume j, k > 1. 
As in Prop. 4.1, we use the decomposition, 



'■ Aj ' h = (^7T'^'^) + 2^ (—>-) =2" A V^(j + l,fc + l,l) + 2- K n-^(j-fc,l), 



which gives 

To begin, we want to show that 



LntJ-l 

lim V n~ l \ip(j + l,k + 1,1)1 = 0. (30) 

3,fc=l 

Proof of (J30j) . Note that 99 = if K = 1, so we may assume if < 1 and if > 1/6. From ([29| . when 
t Ar > 2s and |i — r| > 2s we have 

|p(t, r, 8 )| < 4fi 2 ff (1 - K)s 2 [(t - s) 2H + (r - s) 2H ] A " 2 (t - s) 2H -\r - s) 2 "" 1 

so that 

-2ffif „-i?if- 1 v,HK- 1 



| V (j + l,fc+l,l)|<Cn-^j 
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Recalling the notation Jd from Lemma 3.2.d, we have 

LntJ-l 

n- 1 J2 \<P(j + l,k + l,l)\=n- 1 E Wi + U + MH+n" 1 £ + 1,*+ M)l 

j,k=0 ti,k)£J d V,k)$J d 

'LntJ-1 

< C\nt\n? + C\nt\ 2HK n~^ < Cn~% , 

where we used the fact (which follows from Lemma 3. 2. a and the definition of <p and ip) that \tp(J + 
1, k + 1, 1)| is bounded. 

Hence, (|30|) is proved. It follows from (|30|) that 

[ntj-l 

\<p 3 + 3<p 2 Tp + 3<pi> 2 \ ^ n 0, (31) 

" 3,fc=l 

since <p and V 1 are both bounded. Hence, it is enough to consider 

[ntj-l 

ij(t) = lim J] W - k, If. (32) 



n— >-oo 

3,fc=0 



To evaluate (I52l. we have 



LntJ-1 



8^n 

j,k=0 



— E (li-fc + i|*-2|i-*|* + |j-fc-i|4 



8^n 

i,fc=o 

M-l ., L'"J I ./ -I 



g, I 2 L«*J-i i-i i 3 



8 K n 8 K n 

j — 1 m— 1 



8^31 8 K 71 

j — 1 m — 1 

2 L™*J-! oo 



8 K n 

3=1 rn=3 

where the last term tends to zero since 

[ntj— 1 oo 3 oo 

3=1 "i=J 3=1 

We therefore conclude that T](t) = C^t, where 

1 / °° / v t\ 

1 / _ - /, ^ v 1 1 , , 1 \ I 



C K = 8 + 2 E (fin + 1)* ~ + (m - 1)* 



m= 1 
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This number is approximately 8 ■ □ 

As an immediate consequence of our proof of Theorem 3.1, we have an alternate proof and extension 
of previous results in Gradinaru et al. In [5] , it was proved that ([2J holds for any fractional Brownian 
motion with H > 1/6, that is, the correction term vanishes. Following Remark 3.5, we may conclude 
the following: 

Corollary 4.3. Let B t = {Bf' K ,t > 0} be a bifractional Brownian motion with parameters 1/6 < 
HK < 1. Then on a fixed interval [0, T] and for < s < 1, B satisfies @. 

Proof. Notice that s < 1 implies s 2HK < s^ . With small modifications to the proof of Prop. 4.1, it is 
easy to verify that conditions (i) - (v) are satisfied when HK > 1/6. In accordance with Remark 3.5, 
we want to show that 

[ntj-1 

lim T \/3 n (j,kf\=0. 

n—toc * — ' 
j,k=0 

We may assume K < 1. From Prop. 3.1 of [7], we have that 

E [(B t — B t - S ) 2 ] <Cs 2HK . 

Recalling the notation Jd from Lemma 3.2.d, Cauchy-Schwarz implies for (j, k) € Jd, we have \ f3 n (j, k)\ < 
Cn- 2HK . For (j, k) £ J d , by ((29j) we have 

\<p(j + l,k + 1,1)| < AH 2 K{1 - K) [j 2H + k 2H ] K ~ 2 fH-l k 2H-i 

< c\j - fci" 1 - 2 ^ 1 -^, 

and similar to Prop. 4.1, we have 

m-k,i)\<c\j-k\ 2HK - 2 - 

hence for (j, fc) </ J d we have \/3 n (j, k)\ < Cn- 2HK \j - for 7 = min{l + 2H(1 - K),2HK - 2}. It 
follows that 



[ntj - 1 



< £ |/3„(j,fc) 3 |+ E |/?„(j,fc) 3 | 
0',fc)eJ d (J,k)t£J d 

< J2 Cn- 2HK + Cn- 6ffA ' 2 IJ - fc l 



so |t7(£)| = because i/if > 1/6. □ 
4.2 Extended bifractional Brownian motion 

This process is discussed in a recent paper by Bardina and Es-Sebaiy [1]. The covariance has the same 
formula as standard bBm, but it is 'extended' in the sense that 1 < K < 2, with H restricted to satisfy 

< HK < 1. Within the context of this paper, this allows us to consider values of 1/12 < H < 1/6. 
As in section 4.1, we show computations only for the case HK = 1/6. A result similar to Cor. 4.3 can 
also be shown by modification to the proposition below. 

Proposition 4.4. Let Y = {Y t H ' K , t > 0} be an extended bifractional Brownian motion with parameters 

1 < K < 2, HK = 1/6. Then Y satisfies conditions (i) - (vi), with 9 — 2/3, A = (2H) A ^, and with v, 
7 and r/(t) as given in Prop. J^.l. 
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Proof. Conditions (ii) and (v) are the same as for standard bBm, as shown in Prop. 4.1. In particular, 
the decomposition into <p(t,r,s) and tp(t — r,s) for condition (v) is the same, so it follows that n(t) of 
condition (vi) has the same form. The proofs for conditions (i), (iii) and (iv) require some modifications 
to accept the case K > 1. 

Condition (i). From Prop. 3 of [1] we have 

|E [(Y t - Y t ^ s f} | < s 2HK = si. 

Condition (iii). First, we have 



E [(Y t - Y t ^ s f - (Y t _ s - Y f _ 2s ) 2 ] 

= t 2HK J, [t 2H + (t~ sf H ] K + ^[{t-sf H + {t- 2s) 2 «] K A 

2^ 



< 2t 2HK - ^ [t 2H + (t- s) 2H ] A - 2(t - 2 S ) 2HK + ^ [(t - s) 2H + (t - 2s) 2/r 

f-0 



4ffK / (t + O ai,A_1 -(t-« + 



— -S 



2/. o„\2ffA'-2 



< CV(t- 2s) 

On the other hand 



j^HK 2_ [t 2H + {t _ s?H] K + 2_ [{t s)2H + {t _ 2s)2H] «_2_ {t „ 2sf HK 

> 2(t - s) 2HK - A + 4 2ir] K + 2_ [(; _ 2s )2ff + ( t _ 2s ) 2 "] K - 2 (t - s) 2HK 

= -4FX f (t - s + t^) 2 ^- 1 - (t - 2s + ry) 2 ^- 1 dr? 
Jo 

> -Cs 2 {t~2s) 2HK - 2 1 

hence the term is bounded in absolute value as required, with v = 2 — 2HK = 5/3. 
Condition (iv). 

|E [Y r (F t - y t _ s )] | = i ( [r 2 " + i 2 "] * - [r 2H + (t s) 2H ] K ) + ± (\r - t + s\ 2HK \r t\ 2HK ) 

< _L ( [ r 2H + f 2 H] * _ [r 2H + {t _ sfH] ^ + l_ (Qr _ tl + sf HK |p _ ^ 

We consider two cases for the first term. If t < 2s, then by Fundamental Theorem of Calculus, 

i k 



_L ( [r ™ + t 2 H] K [r2H + {t _ s?H] K^l_^ [r2H + (2sfH] « _ r2HK 

= §<J q V 2H + u] K ~ 1 du<Cs 2H . 

If t > 2s, then 

J_ ( [ r 2H + t 2H] A' _ + {t _ s)2H ^ = J° p ff + {t + u)2H ^ K-l {t + u) 2H-l du 

< CsT 2H ( K -V(t - s) 2 "- 1 < Cs(t - s) 2H -\ 
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In particular, if \r — t\ < 2s then this is bounded by 



Csl " ' r~ s ] c * 2 " 



For the second term, if |j — t\ < 2s, then it easily follows that 

1 

2 K 

and if |r — t\ > 2s, then by Mean Value 



((|r - t\ + s) 2HK - |r - t\ 2HK ) < Cs 2HK < Cs 2H ; 



±_ ((| r _ t \ + s) 2HK -\r- t\ 2HK ) < Cs\r - t\ 2HK ~ l < CsT 2H ^ K -^\r - t\ 2H ~ l < Cs\r - t\ 2H ~\ 



2 K 

In particular, if t < 2s then 



Cs\r - t\ 2 "-' = Cs zn ( y—^ ) < Cs 2H 



Hence, we have shown that 

\E[Y r {Y t -Y t _ s )}\ < 

and so condition (iv) is satisfied by taking A = min {2H, |}, where K e (1,2) implies A > 1/6. □ 



' Cs[{t- s) 2 "- 1 + \r -t] 2 "- 1 ] HT>2s and |r - t\ > 2s 
Cs 2H otherwise 



4.3 Sub-fractional Brownian motion 

Another variant on fBm is the process known as sub- fractional Brownian motion (sfBm). This is a 
centered Gaussian process {Z t ,t > 0}, with covariance defined by: 

R h (s, t)=s h +t h -± [(s + t) h + \s- t\ h ] , (33) 

with real parameter h G (0, 2). Some properties of sfBm are given in [3] and [4]. Note that ft = 1 is 
a standard Brownian motion, and also note the similarity of Rh{t,s) to the covariance of fBm with 
H = ft/2. Indeed, in [4] it is shown that sfBm may be decomposed into an fBm with H — ft/2 and 
another centered Gaussian process. 

Similar to Section 4.1, we discuss only the case ft = 1/3. For ft > 1/3, it can be shown that conditions 
(i)-(vi) are satisfied with 77^) = 0, hence ([2]) holds. 

Proposition 4.5. Let Z = {Z tl t > 0} be a sub-fractional Brownian motion with covariance Q33p and 
parameter ft = 1/3. Then Z satisfies conditions (i) - (vi) of Section 3; hence Theorem 3.1 holds. For 
condition (vi) we have rj(t) = Cut, where 

\ 00 3 

m— 1 

Proof. Condition (i). We have 

E [(Z t - Z t _ s ) 2 ] = R h (t,t)+R h (t-s,t-s) - 2R h (t,t- s) 

= ~t h + \{2t - s) h - ^(t - sf + \{2t - s) h + s h 

= -\ [{2tf - (2t sf] \ [{2t 2s) h (2t sf] + s h . 
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This is bounded in absolute value by Cs h , using the inequality a h — b h < (a — b) h . 
Condition (ii). 

2 h 



|E [Z 2 t - Z?_ s ] 



2 h 

2t h - —t h - 2{t -s) h + —(t- s) h 



By Mean Value this is bounded by 



Cs(t- s) 11 - 1 = Cs(t- s)- 



which implies (ii) with 9 = 2/3. 
Condition (Hi). 

E [(Z t - Z t _ s f - (Z t _ s - Z t _ 2s f] = R h (t, t) - 2R h (t, t-s) + 2R h (t -s,t-2s)- R h (t - 2s, t - 2s) 

= ~^Y~ + (2* - s) h - (2t - is) h + \{2t- As) h 



= -)- [{2tf - 2{2t - sf + (2t - 2s) h ] + \ [{2t - 2sf - 2{2t - 3s)' 1 + (2t - 4s)' 1 ] 



By Mean Value, these terms are bounded in absolute value by 

Cs 2 (2t - 4s)' 1 - 2 < Cs^ +V {t - s)~ v 

for v = 5/3. 
Condition (iv). 

\E[Z r (Z t - Z t _ s )}\ = \R h {r,t) - R h (r,t - s)| 



t h -(t- s) h -±[( r + t) h -(r + t- s) h ] + 1 (|r - t + s\ h - \r - t\ h ) 



du 



Note that the above expression is always bounded by Cs h by the inequality a h — b h < (a — b) h . Hence, 
the bound is satisfied for the cases t < 2s or \t — r\ < 2s. Assuming t > 2s, \i — 1\ > 2s, we have 

t h -(t-s) h -l[(r + t) h -(r + t-s) h ]+±(\r-t + s\ h -\r-t\ h ) 

<h f {t + u) h - 1 du+^ f ( r + t + u) h - 1 du+^ [ (Ir-tl+u) 11 - 1 
J-s 2 J_ s 2 J 

< Cs(t - s)"i + Cs(\r -t\- s)"t. 
For \r — t\ > 2s, we have (|r — t\ — s) > \\r — t\, so 

Cs(t - s)~i + Cs(\r -t\- s)"I < Cs(t - s)*" 1 + Cs\r - t]^ 1 

for A = 1/3. 
Condition (v). 

E [(Z t - Z t _ s )(Z r - Z r _ s )\ = R h (t, r) -R h (t- s, r) - R h (t, r - s) + R h (t - s,r - s) 

= -\ [(t + rf - 2(t + r - s) h + (t + r - 2s) h ] + I [|t - r + s\ h - 2\t - r\ + \t - r - s\ h ] . 
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Assuming that \t — r\ > 2s, by Mean Value this is bounded in absolute value by 

Cs 2 \t-r- s\ h - 2 < Cs 2 \t-r\ h - 2 

since \t — r\ > 2s implies \t — r — s\ > h\t — r\. If ft, < 1, then we take 7 — 2 — h = 5/3, and we have an 
upper bound of 

Cs h+2-h\ t _ r |h-2 = Cs h+i\ t _ r |- 7 . 
Condition (vi). First assume h = 1/3. Referring to condition (v) above, we can decompose fi n {j,k) as 

PnU, k) = ^(j, k, 1) + ^J(j -k,l), 

where u(j,k,l) = -(j + k+2) h +2(j+k+l) h -(j + k) h and^(j-fc,l) = \j-k+l\ h -2\j-k\ h +\j-k-l\ h . 
Note that if>(J — fc, 1) is identical to the ip used in Prop. 4.1, where in this case h = 2HK . Following 
the proof of Prop. 4.2, it is enough to show 

\nt\-l 

hm n- 1 V \u(j,k,l)\=0, (34) 

n— f 00 L — * 

j,k=Q 

so that, similar to ([32]) in the proof of Prop. 4.2, we have 

[ntj-l 
j,fc=0 

where 

= 1 + - ^ ((m + 1)4 - 2mi +(m - 1) 3 J . 

m— 1 

That is, corresponds to the constant Cjc from Prop. 4.2 with K = 1. 

Proof of (f34|) . By Mean Value and the above computation for condition (v), |oj (j, fc, 1)| < C(j + k) 1 
for some 7 > 1. Hence, for each j ' > 2, 

poo 

<C u~<du < C(j-l) 1 -' 1 . 
Jj-i 

It follows that we have 

\nt\ — 1 [nt\-l \nt\-l \ nt\-l 

n- 1 Hj,k,l)\=n- 1 £ (| W (0,M)| + Kl,M)|)+n- 1 £ £ l"0'.M)| 

j',fe=0 fc=0 j=2 fe=0 

[ntj-l 

= n" x £ ([(fc + 2)' l -2(fc + l)' l + fc /l ] + [(fc + 3) /l -2(fc + 2)' l + (fc + l)' 1 ]) 

fe=0 

[ntj-l 

< Cn" 1 +Cn- 1 [nt\ 2 -' f 

which converges to since 7 > 1. □ 
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5 Proof of Some Technical Lemmas 
5.1 Proof of Lemma 3.3 

We may assume t\ = 0. For this proof we use Malliavin calculus to represent AX^ as a Skorohod 

integral. Consider the Hermite polynomial identity x 5 — H${x) + 10H3(x) + 15Hi(x). Using the 
isometry H p (X(h)) = <5 p (/i® p ) (when \\h\\sj = 1) we obtain for each < j < [nt 2 \ - 1, 

AX| = S 5 (df) + 10||AXx||2 2 5 3 (af 3 ) + 15\\AXj_\\i,5(di_). 
With this representation, we can expand 

L«i 2 J-l 

E [/ (S) (^i)/ (S) (^i) AX i A ^i 

j,k=a 



(35) 



into 9 sums of the form 



3,k=0 



(36) 



where p, q take values 1, 3, or 5. By the integral multiplication formula ([5]); and using the Malliavin 
duality ([5]), each term of the form (|36p can be further expanded into terms of the form 



\nt 2 \-l 



C J2 \\^X^\ 5 - p \\AX^\\ 5 - 2 q E f^(X^)f^(X i )S p+q - 



2r ( pf&p— r 



j,k=0 



[nt 2 ]-l 

j,k=0 n 



L 2 ^ 



jjp+q-2r 



(/(5)(Xx)/( 5 )(li)),0f- r ®af 9 - 



where < r < p A q and p,q <E {1, 3, 5}. For < m = p + q — 2r < 10, we have 



D' 



/ (S) (^)/ (5) (^)J = E ^ a (/ (5) (^))^ b (/ (5) (^) 

a-\-b— m 

= E / (5+a) (^)/ (5+fc) (^)eT®e 



a+b— m 

Hence, we expand (|36[) again into terms of the form: 

|nt 2 J-l 

C J2 l|AXi||^ p ||AXi||^*E[/( 5+a )(Xi)/( 5+,, )(Xi 

j,k=0 * 



x ( e 



where a+6=p+g— 2r. With this representation, we are now ready to develop estimates for each 
term. By condition (0), 



E 



f i5+a) {X±)f( 5+h \X h ) 



< 



sup 

v 0<j<|ntaJ 



/( 5 +«)(X X ) 



sup E 

. 0<fc<fn*2j 



/( 5+b )(X A ) 
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and by condition (i), 



sup ||AX X ||^7 P sup \\AXj, < Cn r 



0<j<[nt 2 ] n 0<k<[nt 2 \ 

If a > 1 with a + b = p + q — 2r, by condition (iv) 



e® a ® e r f, 



dT~ r ® d® q - r ) 



< Cn-( a+b - 1 '> x 



ij ,dj 



with a similar term in fc if a = and 6 > 1. Hence, assuming a > 1, each term in the expansion of (|36[) 
has an upper bound of 



Lnt 2 J-l 



c 



i,fc=0 



»" 10-(p+g) 



-(o+fe-l)A 



i±,dj_ 



To show each term has the desired upper bound, first assume r > 1. Then 
by Lemma 3.2.d we have an upper bound of 



» n / a 



< Cn , and 



J,*=0 

which is less than or equal to C\nt2\n~z because A > 1/6. For cases with r = 0, then either a > 1 or 
6 > 1, so without loss of generality assume a > 1. For this case with Lemma 3.2.C we have an upper 
bound of 

LniaJ-l 



J,*=0 



>--(a+b-l)\ 



which is less than C\nt2\n 3 since A < 1/3. 



5.2 Proof of Lemma 3.6 

Without loss of generality, assume a = 0. First we want to show that for each integer < k < b — 1, 

A- 



E 



< C. 



(37) 



Using the Taylor expansion similar to Section 3.2, 

W2+1) - f"{x L ) = (nx^) - f'{X L j) - (f{Xi.) - f"{X L ) 



= f^(X i )AX i + ^(X^AXl + JLfWtf^AXl + B+(j) - B~(j) 



where B+(j),B~(j) have the form C/( 9 >fe-)AA^. Hence we can write, 



E 



E/ (3) (ii)AA, 



< E 



E(mi±i)-r(*x) 



i=0 



-E 
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£/( 5 >(X x )AX| 



2 5 5! 



E 



E J2\ B n(j)\+\ B nU)\ 
3=0 
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We have the following estimates: By condition (0), 



E 



3=0 



< E 



by Lemma 3.3, 



and by Lemma 3.4, 



E 



3=0 



< Cn-s(fc + l)3; 



J2 E \ B nd)\+ E \ B nti)\ < Cn-i{k + lf. 

3=0 

This leaves the AX 3 term. Using the Hermite polynomial identity y 3 = H 3 (y) + 3Hi(y), we can write 



E 



E/ (5) (^l)AX 3 

3=0 



< E 



£/ (5) (*i)<5 3 (df 

3=0 



+ E 



L 2 



For the first term we have 



£./ (5) (^W) 

j=0 



= E e [/^(Xi)/( 5 )(li)^ 3 (af) «5 3 (af) 

fe 3 

= EE 

j,£=0 r=0 



I ;( I E 



fe 3 / \ 2 

Hi!"] E 



-EE 

fe 3 

^ EEE 

j,£=0 r=0 a+b= 
6-2r 



/ (5) (ii)/ (5) (ii)^- 2r (af 3 ~ r ® r af 3 - r ) 

/ (5)(^ i)/ (5) ( ^ ) 



6-2r 



Si Q$ e 



d ± ,de 



For this sum, if r = we use Lemma 3. 2. a and 3.2.b for each pair (a, b) to obtain terms of the form 



E 

3,t=0 







3 




< sup 




sup 




3\t 









< Cn- 1_3A (A; + 1), 



where we use the fact that r = implies a > 3 or b > 3. If r > 1, we use Lemma 3. 2. a and 3.2.d to 
obtain terms of the form 



Cn -(6-2r)A J2 | <Cn-( 6 - 2r ) A -5(fc+l), 



noting that (6 - 2r)A + § > 1. 
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For the other term, we have by Lemma 3.2.b, 

\ 2 



< sup 

0<j<fe 



AX i\\ L , E |e[/ (5) (^)/ (5) (^i) (S 2 ) 



j,t=Q 



<Cn'i E \ D2 f i5) ( x ±)f {5) (X±) 



k 



+ 



E 



/ (5) (li)/ (5) (li) 



dj_,dt 



3,f=0 a+6=2 
fc 



fco / \j=a " / 



< Cn" 1 ^ + I) 1 - 9 + Cn-*(k + if- 26 + Cn- l {k + 1) < C, 
where the estimates follow from Lemma 3.2.C and 3.2.d. Hence, by Cauchy-Schwarz 

k 



E 



3=0 



<C, 



which proves (j37|) . Now we define 



G n (j) = J2fW(X±)AX±, 



k=0 



and by Abel's formula and condition (iii) we have 



E 



b-l 



X AXi 2 / (3) (X i )AX i 

3=0 ^ 



< 



Mi 



6-1 



L 2 



E|G n (6-l)| + ^E|G n (i)l 



3=0 



L 2 



Al 2 



b-l 



< Cn-i. 
5.3 Proof of Lemma 3.10 

Proof of (jM)) . Let = lj and bj- = [ntj\ . By Lemma 2.f .b, 
b,-i 

D 2 Fi = J] D 2 5 3 (/( 3 )(li)9f 3 J 

b,— i 



X {^ 3 (/ (5) (*±)3f) eT + 6 <5 2 (/W(Xi)9| 2 ) 0£ ®£i +6 5 (/< 3 >(Xi)0i) 9f } 
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and 



b k -l 



DF n= E { <53 (/ (4) (^v) 9 f)^+ 3(52 (/ (3) (%)5| 2 )%}. 



(38) 



With these two expansions, it follows that the expectation 



E 



(ui,D 2 Fi®DFl;) s 



consists of terms of the form 

bi-l bj-l b k -l 



p,p—a,i q,q—a,j m,m—ak 



E E E ^ G{p,p')5 r - (g^X^df 1 ) S r > (giiX^A (g 3 (Xm)dp) <T« ( 9i (X L 

r> f — n._- o a 1 =n ■ m rn'—m. L ^ ' ^ 



£V , cV 



8ql_ , 8pl^ 



3— r 2 



\ ri — 1 / \ 3-n 

X (£<*,&£ ) (d*,dp) 

• «•/ £ \ « ™ / \ - — / 

x(di*,dz} \^mL,d^_J^ \pmL,d^_J^ 



r 3 -2 



(39) 



where G(p,p') := f^(Xp_)f^(X^_), r\,T2 take values 1,2 or 3; r 3 ,r^ take values 2 or 3; and each gi 

represents the appropriate derivative of /. Without loss of generality, we will assume that u l n ,D 2 Fl, 
and DF% are all defined over the interval [0, t], and that all sums are over the set {0, . . . , \ nt\ — 1}. Let 
R = n + r 2 + r 3 + n, and note that 6 < R < 12. It follows from Lemma 3. 2. a, 3.2.C, and/or 3.2.d that 



[ntj-l 

E \( i i> 9 $ 

p^p'—O 



si \ d ^)* \ e i>%)» 



\ 3-r 3 / \ r 4 -2 / v 3-r 4 



< E 

[ntj-l 

x E 

p'=0 

< Cn- A , 



ri-l / \3-n / \r 3 -2 / i3-r 3 

£i,9ji) (^d3-,dpj {im,dpj {dm., dp j 



r 2 — l 



d^_,d^_ 



3— r 2 



r-4-2 



3 — T4 



where the exponent A is determined by {n, . . . , r 4 } as follows: First, suppose n = 3. Then by Lemma 
3.2. a and 3.2.C, 



[ntj-l 

E 

p=0 



: -i - <V n ) (^im,dp^ {dm, dp) 



< sup 



r 3 -2 I i 3-1-3 

Em, dp) (dm, dv) 



[nt\ - 1 




E 

p=0 





< (7 n -2A-(r3-2)A-I(3-r3) 

= (7n~ < - ri+r3 ~ 3 ' |A ~^ 6 ~ ri ~ r3 ' ) 
On the other hand, if n = 1 or 2 then by Lemma 3. 2. a and 3.2.d, 

|nt|-l 



E 

p=0 



es_,diL) \da_,di_j {em,dE.J {dm, dp j 



< sup 



sa_, dp 



n-l 



sup 

m,p 



r 3 -2 / \3-r 3 
im,dp ) {dm, dp) 



[nt]-l 

E 

p=0 



d±,dp 



3-n 
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Combining this with a similar computation for the sum over p' , we obtain 

A = X(R — 6) + -(12 - R) = 2 - Q - Aj (R - 6). 

In particular, A = 2 if R — 6 and A = | + A for R = 7. It follows that we want to find bounds for terms 
of the form 



Cn- A 



sup Y, 



3i ,32 ,3a ,34 



(40) 



By repeated use of ([S]), we can expand each product of the form 



into a sum of terms of the form 



Ct-3 



where Cm is a combinatorial constant from ([5]), V&n = Y\t=i 9i{-^-ii)'i eacn a i € {0,1,2}, such that 
A := Yli=i a i — each nonnegative integer bi satisfies 6j < r^; and the exponent M satisfies: 

M = &! + 6 2 + 63 + 64 = -R - 2A 
With this representation, and using the Malliavin duality (|S|), we want to bound terms of the form 



Cn- A sup J2 
p,p' 



E 



D M G(p, P '),* n df^®---®d® bi ) 1 (d^^^Y 1 ••• (d^,dA 



J1>J2,J3,J4 

Consider first the case A = 0. Then M = i? > 6, and each &i = Tj > 1. Hence 



(41) 



Cn- A 



sup ^ 

v,v . . . . 

31,32,33,3- 

< C'n~ A sup 



fl-4 



[ntj— 1 
SUp (££,9j 



4 



By Lemma 3. 2. a and 3.2.c, this is bounded by Cn~ 1 ~ 2X , since A > 1 for all R and R > 6. 

If A > 1, by permutation of indices we may assume that a\ > 1, so (|41[) may be bounded using 
Lemma 3.2.c and 3.2.d: 



[nt]-l 



Cn A sup sup 



E 



x sup 

3,k 



D M G(p, P '),* n d^ 
4-i L«^i 

31,32=0 



< CL^J 3 n _A ~3 sup 

pj 

< LntJ 3 "" 6 , 



1 i\/ 



££ , dj_ ) sup 
I j,k 



,4-1 
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where 



9 = 4+(i?-6 + M)\ - =4+(R-A)(2\-h- 6A. 

o o 

Since A < R/2, R > 6, and A > 1/6, we have > 3 for all cases except when R — 6, A — 3. This case 
has the form, 



Cn~ 2 sup J2 \ E [ G (P>P')*n}(dn,diz) ■■■(di3,d l 

p p' ....I \ n n / \ji7 

JU2J3J4 



< Cn 2 sup 



31,32 



33,34 



ai — X i \ct2 I \aa 



\nt\ - 1 



3l,h=0 



E \\ d if d i 



< C[nt\n~ 



\_nt\-l 
sup ^ I (d± , 9*; 



sup 

3,k 



[ntj-l 

33 ,34=0 



by Lemma 3.2.d. 

Proof of ([2~5j) . For this term, we see that 
consists of terms with the form 

bi-l bj-l b k -l bi-l 

E E E E : 

p,p'=a,i ji,32=aj 33,34=0* js,3e= a e 



where each rj € {2, 3} and G(p,p'),gi(x) are as defined above. As with ([24)) above, we assume that all 
components are defined over the time interval [0,f] for some t < T. As above, let i? = X};=i r i' an d 
note that for this case 12 < i? < 18. Similar to the above case, we obtain 



E 



ri — 2 / y3-ri 



re — 2 / y 3-re 



<Cn- A ', 



where A' = 2 - (§ - A)(i? - 12). It follows that, similar to (gUJ), we want to obtain bounds for terms of 
the form 



[nt] - 1 

Cn- A ' sup jr 

PlP ' ji,...,j6=0 



E 



i=i 



Using ([6]) and the Malliavin duality as before, we obtain terms of the form 



[ntj-l 

CrT K sup 



p,p' 



ji,...,j 6 =0 



D M G(p, P ),S M l^ndf* 1 ® • • • ® 3 



n |(%'% 

{jlj — iJ's} 



(42) 

where *$> n — Yli=i 9i(Xj±), each on and each 6i take values from {0,1,2,3}; and the product includes 
all 15 possible pairs from the set {ji, . . . such that A := X)i=i a i — -R/2- As in the above case, for 
each R we have M and A satisfying M = X)i=i anc ^ M = R — 2 A. 
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In the product 



D M G(p,p),S M U n dfi®..-®df°)\ TT (duA 



{ii,...,j 6 } 



(43) 



each of the indices {ji, . . . must appear at least once. Note that by Lemma 3. 2. a we have (possibly 
up to a fixed constant) 



sup 


fa' 9 *)* 


< sup 




0<j,k<[nt\ 




0<j,p< [rat J 





and by Lemma 3.2.C and 3.2.d we have 

[ntj-l \nt\-l 

sup 2_, \\9±,dh) < sup 2_, 

0<k<|nt| j=0 1 X " 0<P<|ntJ i=0 ' " 

Hence, we may conclude that ((32]) contains terms less than or equal to 

(%>dik) ^££,9j3.)^ (^f'^ia.)^ ' 
and, by Lemma 3.2.c, (|4"2")l is bounded in absolute value by 



M-3 



sup 



dj , dk 



" I S3 



[ratj-l 
SUp ^ |^££,(9i 
" 1 3=0 



[ntj-1 

Cn~ A ^ sup (es_,d± 

3i,h,je=0 PJ 

< C[nt\ 3 n- e ' ', 

where, using the fact that ii = M + 2A, 

1 ^4 1 

6' = 2- (- - X)(R - 12) + (M-3)A+ - = 6+ (R- A)(2X - -) - 15A. 

Observe that 0' > 3 whenever R — A > 6. The case i? — A = 6 occurs only when R = 12, A = 6, and 
M = 0; so in this case we have an upper bound of 



[ntj-l 
CrT 2 2_, SU P 



7 k 

34,35,36=0 J ' 



dj_,dk 



\nt\-l 

sup ^ 



J=0 

< C[nt\ 3 n- 2 -% < CrT 1 . 
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